Excited state dynamics in DNA double helices 
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Recent ultrafast experiments have implicated intrachain base-stacking rather than base-pairing 
as the crucial factor in determining the fate and transport of photoexcited species in DNA chains. 
An important issue that has emerged concerns whether or not a Frenkel excitons is sufficient one 
needs charge-transfer states to fully account for the dynamics. In we present an SU{2) ® SU{2) 
lattice model which incorporates both intrachain and interchain electronic interactions to study 
the quantum mechanical evolution of an initial excitonic state placed on either the adenosine or 
thymidine side of a model B DNA poly(dA).poly(dT) duplex. Our calculations indicate that over 
several hundred femtoseconds, the adenosine exciton remains a cohesive excitonic wave packet on 
the adenosine side of the chain where as the thymidine exciton rapidly decomposes into mobile 
electron/hole pairs along the thymidine side of the chain. In both cases, the very little transfer to 
the other chain is seen over the time-scale of our calculations. We attribute the difference in these 
dynamics to the roughly 4:1 ratio of hole vs. electron mobility along the thymidine chain. 



Give the importance of DNA in biological system and 
its emerging role as a scaffold and conduit for electronic 
transport in molecular electronic devices, DNA in 
its many forms is a well studied and well characterized 
system. What remains poorly understood, however, is 
the role that base-pairing and base-stacking plays in the 
transport and migration of the initial excitation along the 
double helix. 0,13 Such factors are important since the 
UV absorption of DNA largely represents the weighted 
sum of the absorption spectra of it constituent bases 
whereas the distribution of lesions formed as the result of 
photoexcitation are generally not uniformly distributed 
along the chain itself and depend strongly upon sequence, 
suggesting some degree of coupling between bases. Q 

Recent work by various groups has underscored the 
different roles that base-stacking and base-pairing play 
in mediating the fate of an electronic excitation in DNA. 
0,13 Over 40 years ago, Lowdin discussed proton tunnel- 
ing between bases as a excited state deactivation mech- 
anism in DNA 4] and evidence of this was recently re- 
ported by Schultz et al. 5] In contrast, however, ultrafast 
fluorescence of double helix poly (dA). poly (dT) oligomers 
by Crespo- Hernandez et al.,2| and by Markovitsi et al. 
[2| give compelling evidence that base-stacking rather 
than base-pairing largely determines the fate of an ex- 
cited state in DNA chains composed of A and T bases 
with long-lived intrastrand states forming when ever A is 
stacked with itself or with T. However, there is consider- 
able debate regarding whether or not the dynamics can 
be explained via purely Frenkel exciton models 0, Q, 
or whether charge-transfer states play an intermediate 
role, 

Here we report on a series of quantum dynamical calcu- 
lations that explore the fate of a localized exciton placed 
on either the A side or T side of the B DNA duplex 
poly(dA)io.poly(dT)io. Our theoretical model is based 
upon a SU{2) (g) SU{2) lattice model consisting of lo- 
calized hopping interactions for electrons and holes be- 




FIG. 1; Three dimensional structure of poly(dA)poly(dT) in 
the B DNA double helical form. This is considered to be the 
average structure in water. Bottom: Equivalent lattice model 
showing the connectivity and associated transfer terms. 



tween adjacent base pairs along each strand (taj) as well 
as cross-strand terms linking paired bases (hi) and "di- 
agonal" terms which account for the tt stacking inter- 
action between base j on one chain and base j ± 1 on 
the other chain (r^) in which rj denotes coupling in 
the 5 '-5' direction and coupling in the 3 '-3' direc- 
tion. Fig. n shows the three-dimensional structure of 
poly (dA)io. poly (dT) 10 and the topology of the equiva- 
lent lattice model. Taking link in this figure as a specific 
electron or hole hopping term, we arrive at the following 
single particle Hamiltonian, 

hi^Y^ CjV'jV'j + tjii'l+i^j + i'li'j+i)) + hjijj^j 
j 
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where i/ij and '^j are SU{2) spinors that act on the 
ground-state to create and remove an electron (or hole) 
on the jih adenosine or thymidine base along the chain. 
The 7 operators are the 2x2 Pauli spin matrices with 
V'j = 7iV'] S'Ud. 7+ + 7_ — 7i providing the mixing be- 
tween the two chains. Taking the chain to homogeneous 
and infinite in extent, one can easily determine the en- 
ergy spectrum of the valence and conduction bands by 
diagonalizing 

~ / e, + 2i,cos(g) h + r+ e-''" + e+^'i \ 
^ \h + r+e+''i + r-e-''i eb + 2ti, cos{q) > 

where ea.b and ta,b are the valence band or conduction 
band site energies and intra-strand hopping integrals. 
When = rj, Eq. [21 is identical to the Hamiltonian 
used by Creutz and Horvath ^3 to describe chiral sym- 
metry in quantum chromodynamics in which the terms 
proportional to r are introduced to make the "doublers" 
at g oc TT heavier than the states at g cx 0. 

The coupling between the conduction and valence 
bands is accomplished by introducing short-ranged 
Coulomb and exchange interactions as well as "dipole- 
dipole" terms which couple geminate electron-hole pairs 
on different sites. 

H{12) = hi + h2 + J2 Vm,r.AlAn (3) 

m,n 

where the are spin-symmetrized composite opera- 
tors that create or remove singlet or triplet electron/hole 
pairs in configuration |m) = \iejh) where Vmn = 
-(mf.nh\\nf.mh) + 2S.q n(rn^nh \\n^rnh) where S" = 1,0 is 
the total spin. jilElElEi 

The single particle parameters are taken from Anatram 
and Mehrez as determined by computing the Coulomb 
integrals between HOMO and LUMO levels on adja- 
cent base pairs with in a double-strand B DNA sequence 
using density functional theory (B3LYP/6-31G) [isj . 
For (dA).(dT), the intrachain electron and hole transfer 
terms are -0.023 eV and -0.098 eV for the thymidine chain 
and -1-0.024 and -)-0.021 eV along the adenosine chain. 
The interchain terms are hj = 0.063 eV, = — 0.012eV, 
and rJ = — 0.016eV for the electron and hj = 0.002eV, 
r+ = -O.OOTeV, and rJ = O.OSOcV for the hole with 
site energies Se — — 0.931eV and Sh — — 6.298eV for the 
thymidine chain and Se — 0.259eV and Eh = — 5.45eV 
for the adenosine chain. It is important to note that the 
asymmetry introduced with rJ gives directionality 
between the 3'- and 5'- ends of the chain. 

For simplicity, we take the on-site Coulomb interaction 
J = —2.5 cV and the on-site exchange interaction to be 
K ~ l.OeV for both the purines and pyrimidines. We 
assume these interactions to be local since the distance 
at which the Coulomb energy between an electron/hole 
pair equals the thermal energy in aqueous ionic media at 
300K is on the order of the base-stacking distance. These 



we set as adjustable parameters to tune the predicted ab- 
sorption spectrum of our model to reproduce the experi- 
mental UV absorption spectra. Lastly, we estimated 
the coupling between geminate electron/hole pairs on dif- 
ferent bases {nenhWrn^mh) via a point-dipole approxima- 
tion by mapping the tt — tt* transition moments onto the 
corresponding base in the B DNA chain. As in the 
electron and hole hopping terms, the local dipole-dipole 
couphng terms read: = —0.099 eV between adjacent 
base pairs, = 0.0698eV along the adenosine chain, 
dj[ = 0.143eV along the thymidine, dr = -0.006eV 
for the 5 '-5' diagonal coupling and (i+ = — 0.013eV for 
the 3'-3' diagonal coupling. The use of the point-dipole 
approximation in this case is justified mostly for conve- 
nience and given the close proximity of the bases, multi- 
pole terms should be included in a more complete model. 
|l7| As a result, the matrix elements used herein provide 
an upper limit (in magnitude) of the couplings between 
geminate electron-hole pairs. Most importantly, however, 
the point-dipole approximation provides a robust means 
of incorporating the geometric arrangement of the bases 
into our model. Having briefly discussed the model 
and how we determined the parameters, we move onto 
discuss the result of our dynamical calculations. 

We consider the fate of an initial singlet electron/hole 
pair placed either in the middle of the thymidine side 
of the chain or the adenosine side of the chain (i.e. an 
exciton). We assume that such a configuration is the 
result of an photoexcitation at the appropriate photon 
energy (4.87 eV for the thymidine exciton and 5.21 cV 
for the adenosine exciton respecitively) and based upon 
the observation that the UV absorption spectra largely 
represents the weighted sum of the UV spectra of the 
constituent bases, such localized initial states are Since 
these are not stationary states, they evolve according to 
the time-dependent Schrodinger equation. Over a short 
time-step, bt this is easily computed using the Tcheby- 
chev expansion of the time-evolution operator |l8| 

*(i + (5t) = exp[-iH(12)(5t/fi,] *(i) 
_ 

= e"'^'*/''^a„(a)r„(-*#)|*(i) (4) 

ts.E — Emax — Emin- whcrc Tn{x) are Chebychev poly- 
nomials. 

r„+i(x) = 2xTn{x) - Tn-i{x) (5) 

with Tq{x) — 1 and Ti(x) = x. H is Hamiltonian oper- 
ator shifted and scaled so that its eigenvalues lie within 
E £ [—1,1], and J„(a) the nth spherical Bessel function 
with a = AEt/2T%. The advantage of this approach is 
that converges rapidly, preserves norm, and is computa- 
tionally efficient since it involves at most matrix-vector 
multiplications . 

In the top two frames of Fig. [21 we show the transient 
probability for finding an exciton placed on the adenosine 
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FIG. 2: Top; Time evolution of the exciton density for an 
initial excitation placed on the adenosine (left) and thymidine 
(right) chains). Bottom: The corresponding net charge on 
a given site following excitation of the adenosine (left) and 
thymidine (right) chains. 

(left) or thymidine (right) chain at time t — in some 
other excitonic configuration along either the adenosine 
or thymidine chain at some time t later. In both cases, 
negligible exciton density is transferred between chains 
and the excitons rapidly become delocalized and scatter 
ballistically down the DNA chain. 

There are some striking differences, however, between 
the exciton dynamics in adenosine versus those in thymi- 
dine. First, in comparing the matrix elements, one 
easily concludes that the exciton mobility along the 
thymidine chain is considerably greater than the mobil- 
ity along the adenosine chain. This can be see in Fig 2. 
comparing the time required for an excitonic wavepacket 
to reach the end of either chain. In adenosine, the exciton 
travels 5 base pairs in about 25 fs where as an exciton 
along the thymidine chain covers the same distance in 
about 10 fs. This factor of two difference in the exci- 
ton velocity is commensurate with the «1:2 ratio of the 
: cf^ intrachain excitonic couplings. 

Secondly, we note that the adenosine exciton remains 
qualitatively more "cohesive" than the thymidine exciton 
showing a number of ballistic traverses up and down the 
adenosine chain over the 300 fs we performed the calcu- 
lation. One can also note that the exciton velocity in the 
5 '-3' direction is slightly greater than in the 3-5' direction 
as evidenced by the exciton rebounds off site Ai slightly 
sooner than it rebounds from site Aiq. This is due to the 
asymmetry introduced in by the and terms. All in 
all, one can clearly note a series of strong recurrences for 
finding the adenosine exciton on the original site ^5 ev- 



ery 100 fs. The thymidine exciton dynamics are far more 
complex as the exciton rapidly breaks apart. While few 
recursions can be noted, however, after the first ballistic 
traverse, the thymidine exciton no longer exists as a cohe- 
sive wavepacket and is more or less uniformly distributed 
along the thymidine side of the chain. 

The excitonic dynamics only tell part of story. In the 
lower two frames of Fig. |21we show the net charge taken 
as the difference between the hole density and electron 
density on a given base. Note the difference in scale be- 
tween the bottom left and bottom right figures. In the 
case where the initial exciton is on the adenosine chain, 
very little charge-separation occurs over the time scale of 
our calculation. On the other hand, when the exciton is 
placed on the thymidine chain, the exciton almost imme- 
diately evolves into a linear combination of excitonic and 
charge-separated configurations. What is also striking is 
that in neither case do either the electron or hole transfer 
over to the other chain even though energetically charge- 
separated states with the electron on the thymidine and 
the hole on the adenosine sides of the chain are the low- 
est energy states of our model. 0| It is possible, that 
by including dissipation or decoherence into our dynam- 
ics, such relaxation will occur, however, on a time scale 
dictated by cross-chain transfer terms. For the coupling 
terms at hand, electron or hole transfer across base pairs 
occurs on the time scale of 3-4 ps. 

The difference between the excitonic dynamics follow- 
ing excitation of A vs. T can be quantitatively noted by 
comparing the curves shown in Fig. |31 where we compare 
the projection of the time-evolved state onto the exci- 
tonic configurations of the chain on which the exciton was 
placed {Paa and Ptt) compared to the projection onto 
the excitonic configurations of the other chain {Pat and 
Pta)- For the case in which the adenosine chain was ex- 
cited, approximately 75% of the total probability density 
remains as excitonic configurations along the adenosine 
chain. In stark contrast, only about 40% of the initial 
thymidine exciton remains excitonic along the thymidine 
side of the chain. The reason for the remarkable differ- 
ence between the two chains stems from the difference 
in electron and hole mobility along the thymidine chain. 
Indeed, comparing the electron and hole hopping terms 
given above, th/t^ ~ 4 for the thymidine chain compared 
to th/te ~ 1 for along the adenosine chain. This is mani- 
fest in the lower right panel of Fig. 2 where we see almost 
immediately a negative charge remaining for a few fs on 
the site where the initial excitation was placed. 

The results described herein paint a similar picture to 
that described by recent ultrafast spectroscopic investiga- 
tions of (dA).(dT) oligomers in that the initial excitonic 
dynamics is dominated by base-stacking type interactions 
rather than by inter-base couplings. Interchain transfer is 
multiple orders of magnitude slower than the intrachain 
transport of both geminate electron/hole pairs as exci- 
tons and independent charge-separated species. Indeed, 
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for an exciton placed on the adenosine chain, our model 
predicts that exciton remains as a largely cohesive and 
geminate electron/hole pair wave function as it scatters 
along the adenosine side of the chain. Our model also 
highlights how the difference between the mobilities in 
the conduction and valence bands localized along each 
chain impact the excitonic dynamics by facilitating the 
break up of the thymidine exciton into separate mobile 
charge-carriers. In the actual physical system, the mo- 
bility of the free electron and hole along the chain will 
certainly be dressed by the polarization of the medium 
and reorganization of the lattice such that the coherent 
transport depicted here will be replaced by incoherent 
hopping between bases. 

In conclusion, we present herein a rather compelling 
model for the short-time dynamics of the excited states in 
DNA chains that incorporates both charge-transfer and 
excitonic transfer. It is certainly not a complete model 
and parametric refinements are warranted before quanti- 
tative predictions can be established. For certain, there 
are various potentially important contributions we have 
left out: disorder in the system, the fluctuations and vi- 
brations of the lattice, polarization of the media, dissipa- 
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tion, decoherence, etc. These we recognize as lacunae in 
our model. 
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FIG. 3: Net probability, Pif for the initial exciton to remain 
an exciton on the initial chain [Paa&^Ptt) or to be transferred 
as an exciton the other chain [Pat&c-Pta) 
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